Comment on "Dynamic properties in a family of competitive growing models" 
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The article [Phys. Rev. E 73, 031111 (2006)] by Horowitz and Albano reports on simulations 
of competitive surface-growth models RD+X that combine random deposition (RD) with another 
deposition X that occurs with probability p. The claim is made that at saturation the surface width 
w{p) obeys a power-law scaling w{p) oc 1/p^ , where 5 is only either 5 = 1/2 or 5 = 1, which is 
illustrated by the models where X is ballistic deposition and where X is RD with surface relaxation. 
Another claim is that in the limit p — > 0^, for any lattice size L, the time evolution of w{t) generally 
obeys the scaling w{j), t) oc {L" /p^)F{p'^^t/ L''), where F is Family- Vicsek universal scaling function. 
We show that these claims are incorrect. 

PACS numbers: 05.40.-a,68.35.Ct 



In Ref.fl]] the following scaling ansatz is proposed: 



(1) 



where w{p, t) are time evolutions of surface width in com- 
petitive growth models RD+X when a random deposition 
(RD) process is combined with process X, and p e (0; 1] 
is the selection probability of process X. The function 
F{-) represents Family- Vicsek universal scaling. The an- 
szatz dH) has been studied previously ^, ^Sj, [j| by ex- 
amples where X represented either Kardar-Parisi- Zhang 
or Edwards- Wilkinson universal process. The new claim 
that is being made in Rcf. 1] is that a nonuniversal and 
model- dependent exponent 5 in Eq.(IT]) must be only of 
two values, either (5 = 1 or ^ = 1/2, for models stud- 
ied in Ref.[l|. To show that this claim is not correct 
we performed (1 -I- 1) dimensional simulations of RD-fX 
models when X is ballistic deposition (BD) and when X is 
random deposition with surface relaxation (RDSR) , and 
performed scaling in accordance to Ref. l]. Our results 
are presented in figs[l][31 

Our data have been obtained on L site lattices [L is 
indicated in the figures) with periodic condition, start- 
ing from initially flat substrates, and averaged over 400 
to 600 independent configurations. The time t is mea- 
sured in terms of the deposited monolayers. Simulations 
have been carried up to t = 10^, and the surface width 
at saturation has been averaged over the last 5000 time 
steps. The data sets are for ten equally spaced selec- 
tion probabilities p from p = 0.1top=l, where p = Q 
would be for RD process with no X present, and p = 1 
is for process X in the absence of RD. The data have 



been scaled in L with the theoretical values of univer- 
sal roughness exponent a and dynamic exponent z of 
the universality class of process X. The RDSR algorithm 
used in our simulations is given in Ref.|5(| (Sec.5.1). The 
BD algorithm used as X=BD1 is the nearest- neighbor 
(NN) sticking rule found in Ref.'sl (Sec.2.2), and the BD 
algorithm used as X=BD2 is the ncxt-nearest-neighbor 
(NNN) sticking rule found in Ref.g (Sec.8.1). 

Saturation. Saturation data (fig[I|) show that in spe- 
cial cases an approximate power law w{p) oc 1/p^ may 
be observed. However, this is not a principle. Even if 
the data can be fit to the power law in p only one of 
our examples shows a reasonable fit with ^ w 1 (seen in 
figdJi). When X=BD1 the data in figdJa show 5 < 1/2. 
The other two examples shown in figUdefy a linear fit. 
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FIG. 1: (color on line) Interface width at saturation in the 
RD-I-X model vs the selection probability p of process X. (a) 
X is RDSR: the case when both RD and RDSR deposits are of 
unit height (diamonds, RDSRl; L — 500); and, the case when 
RDSR deposits are of unit height and RD deposits are of twice 
that height (squares, RDSR2; L = 100). (b) X is BD: the case 
of the NNN rule (circles, BD2); and, the case of the NN rule 
(triangles, BDl). In RD-I-BD simulations L = 500. Solid line 
segments connecting data points (symbols) are guides for the 
eye. The dashed lines give reference slopes. 
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FIG. 2; (color on line) Time-evolutions w^{p,t) in RD+BDl. 
(a) Scaling in p after Ref. jl]. The insert shows the scaled 
initial transients, (b) Evolution curves before scaling. The 
dashed line is the RD evolution for p — 0. In all models when 
the simulations start from a flat substrate w{t) must pass an 
initial transient before universal scaling can be measured. The 
initial transients in part (b) follow RD universal evolution. 
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FIG. 3: (color on line) Time-evolutions w^{p,t) in RD+BD2. 
(a) Scaling in p after Ref.[l||. The outcome of this scaling 
is summarized in the insert, (b) Data before scaling. The 
dashed line is the RD evolution for p — 0. 

In these cases there is no power law of the type claimed 
in Ref. Jj . This absence of power-law scaling in p is also 
evident in fig. 4 of Ref.[l|. 

The RD limit. Another claim of Ref.[l| is that Eq.((lD 
with the power-law prefactors (where (5 = 1 or 1/2) 
would prevail in the RD limit of p ^ 0+ , and that such 



a scaling would be universal. We tested these claims 
in simulations of RD-I-BD models and found the evi- 
dence to the contrary (figsl^lEl)- order to prove the 
absence of power- law scaling via Eq.([T|) in the RD limit 
we present in figsl2)3-[3j3 the original w'^{p,t) data be- 
fore scaling. These original data show that parameter p, 
p € (0; 1], assigns an order in the set of all curves w'^{p, t) 
in such a way that u;^(l, t) is the lowest lying curve, and 
at p = the initial transients become the RD univer- 
sal evolution ?«|j^(0, i) oc t. The region between the 
boundaries w'^{l,t) and wj^^{0,t) is densely covered by 
the curves 'w'^{p, t) because p takes on continuous values. 
The pattern shown in figsl^b-Eb for p G [0.1; 1] extends 
down to values that are infinitesimally close to p = 0, i.e., 
to the entire range of p. If the simulations are stopped 
at infinitesimally small p' the width w'^{p',t) is always 
the highest lying curve in figsl^b-Eb. In other words, the 
smaller the p' the higher the saturation value of w'^ip' , t). 
But there is no bounding highest curve w'^{p',t) in this 
set since the boundary w^{0, t) is the RD evolution. This 
order is reversed under the scaling of Eq. ([1]) when we set 
,5 = 1/2, following Ref.Q. The outcome of this scaling is 
seen in figsI2^-[5^: the boundary w^(l, t), i.e., the lowest- 
lying curve in figsl^b-lSb, is mapped onto the highest-lying 
curve in the image of this scaling seen in figsl2^-[3^; and, 
a higher- lying curve w'^{p, t) before scaling in figsl2)3-[5)3 is 
mapped onto a lower-lying curve after scaling in figs|2^- 
[5^. In this scaling, the initial transients become ever 
longer as p becomes ever smaller and closer to p = 0, as 
seen in the insert in figH^. For any range of p, also in 
the limit p — > 0+ , the image of this scaling demonstrates 
no data collapse. This image is shown in the insert in 
figlSK. Hence, for RD-hBD models Eq.^ with (5=1/2 
does not produce data collapse. 

In some instances of model X, however, Eq.([T]) can 
give an approximate data collapse [3, 0] but then i5 is not 
restricted to the two values postulated in Ref.Q. For 
example, for the RD-fBDl model such scaling can be 
obtained with S « 0.41 3] (note, 0.4 < S < 0.5 is seen 
in fig[T]3). But for the RD-fBD2 model there is no value 
of S that would produce data collapse when nonuniversal 
prefactors in Eq.((T]) are expressed as a power law p^ . We 
have demonstrated that such scaling does not generally 
exist and if occasionally it is observed it is a property of 
particular model. 

In summary, the form of the nonuniversal prefactors 
as seen in universal Eq.([T]) is a fit and is not a principle. 
The exponent 5 in Eq.([T]) is model dependent, and the 
prefactor that enters may have other forms than p^ . 
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